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SUHKftlg 

The  foUoving  report  discusses  ttie  use  of  similaritgr  flcsrs  in 
certain  ^f^rodynamic  Impact  problems,  as  a  basis  for  app*oxiM,te  aMlytlcal 
solutions.  M'ter  a  qualitative  discussion  of  various  ways  ;  ■which  s\^h 
solutions  may  be  obtained,  a  single  derl-vation  is  obtained  ?  a  class  of  one- 
para^ter  similarity  flows  for  problems  in  1,  2  and  3  dii^nsions,  >&leh 
includes  sidjstantially  all  those  used  previously,  and  aHoro  tl»  superposition 
of  eonst«it  uniform  actions.  'Hie  perturl»tion  eqiaitions  relate  the 

asymptotic  fom  of  tl^  actual  solution  to  stab  similarity  flows  (under 
scsRe  conditions)  are  explicitly  formulated  and  the  equations  of  ■^e 
chaiacberlstic  surfaces  which  determlrfe  their  bel^vlor  are  derived. 
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I.  IM1K)DUCTI0K 

fecause  of  their  non-linear  form,  the  eqiiatlons  of  hydrodynaiaies 
(aM  the  pisDhlems  they  represent)  are  not  auienaMe  to  gemratl  analytical 
solutions.  the  other  hand,  are  often  so  cpaplex  that  direct 
nmerieal  solutions,  if  at  all  possible,  beccsme  uneconomic  and/or  un¬ 
convincing  vi-tiiout  SOBM  kiM  of  analytical  guidance.  Bils  is  particularly 
true  in  •^e  case  of  problems  involving  shocks,  where  the  singular  nature 
of  -ttie  flow  may  lead  to  serious  nisnerlcal  difficulties  if  its  implications 
for  the  solution  are  not  properly  uMer stood. 

It  therefore  becomes  necessary  to  utilise  those  analytical  methods 
that  are  available  to  provide  at  least  a  qualitative  understanding,  if 
possible,  of  the  solution  of  the  problem  posed,  and  even  if  not,  a  guide 
to  Ihe  niar^ical  approach.  Ihere  are  a  variety  of  st^h  techniques 
available  in  ths  various  brenches  of  ^drodynamics,  and  they  are 

by  no  means  always  aideqvtate,  th^  often  provide  the  first  step  of  an 
approximate  analytical  appsoach  ^ich  meets  the  need  described  above. 

These  methods  involve  one  or  more  of  the  following; 

1.  the  use  of  special,  solutions,  e.g.,  e<»iformal 
mapping  in  stationary  iiKOmiaressible  flows, 
legeraire’s  transformation  in  more  ^neral  problems, 
and  ’’similarity"  (bonsslogical)  solutions  in  p^blems 
involving  explosions  and/or  boundary  layers; 

2.  perturbation  techniques  around  such  solutions  (in 
particular  singular  perturbation  methods  for 
boundsuy  layer  p*oblems); 'and 

3.  integral  methods  based  on  ttie  conservattvc  properties 
of  the  flow.  (These  seem  partieiilarly  applicable  to 
problems  involving  shocks.) 
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I  In  tile  case  of  hyperveloeity  itapaets,  one  is  confronted  witii 

I  F'oblems  which  lack  cOTiplete  syinmetiy  and  involve  severe  discontinuities 

I  (shocks),  as  well  as  regions  of  sacoth  f]k>w,  and  so  we  may  expect  a  need 

I  for  all  tile  met^sds  cited  above  in  order  to  provide  a  satisfactory 

I  analytical  fouadation  for  understanding  both  tiie  experimental  and 

I  numerical  results.  Despite  tiieir  general  difficulty,  Imiact  probleos 

I  have  certain  characteristics  i*lch  suggest  that  tiie  above  cited  methods 

I  are  particularly  applicable.  The  initial  and  boundary  coMitions  ai-e 

I  usually  sufficiently  localized  in  si»ce  and  time  as  to  effectively 

I  a^roxiaate  an  explosl<»,  and  a  correspondingiy  similar  "similfia^ity” 

I  solution.  ISie  deviation  of  tiiese  coalitions  from  those  necessary  for 

I  a  similarity  solution,  as  well  as  the  stroag  deviatiec  near  the  * 

I  shock,  is  confined  to  a  relatively  small  region  of  space  and  time,  making 

I  possible  effective  perturbation  methods  (albeit  of  a  singular  nature  in 

I  the  second  case). 

I  Ihere  are  thus  several  approaches  available  for  jroviding  an 

I  approximate  analytical  trammTk  for  uMerstaiidlng  and  furtoer  developing 

I  quantitative  (i.e.  generally  nisaerical)  solutions  of  hyperveloeity  imptet 

j  tybblems.  The  first  step  in  each  of  tiiese  is  the  est^li  steaent  of  so^ 

I  suitable  "approxisating”  similarity  solution.  This  may  then  be  used  In 

%  various  vsys t 

I  1.  as  the  first  stage  in  an  approximation  or  pertui^tion 

I  ]p*ocedure  in  some  parameter  of  the  problem  (e.g.  the 

I  shock  strength),  (Sakural)|^5^ 

1  2.  as  a  general  form  (with  variable  parameters)  in  an  Integral 

I  method,  in  which  the  paroaeters  are  determined  phenOTeno- 

j  log:'cally  fr<^  the  conservation  laws  (Chemyi)i^^^ 

f  (see  also  5telchMtnn^^^  and  Stiart^^^) 

I  3*  as  an  asymptotic  solution,  in  ^Ich  case  it  must  be 

I  shown  tiiat  the  initial  aM  boundary  variations  decty 

I  suitably  as  the  process  develops* 
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(l)  and  (2)  above  have  received  a  good  deal  of  attenticai  in  the 
cases  of  point  explosions  and  hypersonic  flow,  in  laost  cases  wiib  a 
certain  amount  of  success.  In  tte  case  of  hypervelocity  impact  the 
probleais  are  more  difficult  (largely  because  of  eq^uation-of -state 
Hmitations),  nevertheless  some  progress  has  been  mde;  method  (l) 
has  been  applied  by  Rae)^^  using  results  of  Sakurai^^^  and  Oshima  in  the 
explosive  case,  and  a  version  of  method  (2)  has  been  used  by  Raizer.^^^ 

Ihe  third  me'^od  Indicated  above  does  not  seem  to  have  been  used#  It 
involves  ra-Uier  more  nathecatlcal  profundity  even  in  its  formulation,  but 
if  s^KJcessfuUy  carried  out  iffovldes  a  more  satisfactory  conceptual 
framework  for  both  of  the  other  uathods  cited,  as  well  as  for  new  caies 
which  may  be  developed.  It  is  not  proposed  to  carry  such  a  rather 
difficult  program  through  here,  ^hat  will  be  done  below  is  to  give  a 
unified  derivation  of  the  various  one  dimensional  siailarily  solutiorm 
used  by  Raizer  aM  others,  and  briefly  to  formulate  the  mathematical 
probl<^  of  the  a^mptotlc  behavior  of  actual  solutions. 

He  effective  utilizatifSi  of  a  similarity  flow  (or,  indeed,  ary 
©"tiler  flow)  as  a  first  step  in  an  approximation  procedure  depei^s  strongly 
on  its  stmplici-^  and  perspicuity.  In  all  cases  of  interest,  both  for 
explosions  and  impacts,  at'tention  has  been  restricted  to  one -dlmenslo lal 
(or  more  strictly  speaking,  one -parameter)  solutions,  -though  the  ^sical 
problem  itself  nay  of  course  be  multi-dimensional.  Hiis  restrictiai 
Implies  a  subs-fcantlal  degree  of  geometric  symmetry,  more  specifically, 
limitation  to  plane,  s^erically  symmetric,  and  axial  or  crans^’■erse 
cyliMrically  symmetric  problems.  Ihe  question  of  any  deviations  from 
such  symsetiy  will  not  even  be  touched  on  hero. 
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II.  SagLE  SIMILARITY  FLOJS 

Che  thus  considers  the  basic  hydrodynamic  equations  (fo: 
gas)  in  the  fora 

gf  +  ^  •  pv  ^  0 


r  a  perfect 


dv 


^  3t  ^  (v  .  7v)  «  -  7p 


iffi 

¥, 


dt 


(pp*^)  +  V  •  7(pp"^)  = 


0 


Fouowlng  .he  genem  method  of  MiehalC?)  Morgen'®)  end 

>■,  me  seerohes  for  e  one  peraneter  set  of  transformetions 

t  -  bt 


Z  -* 
r  -•  b  r 


Z  k8  "* 

V  b  V 

P  -  p 

p  -♦  b^p 


irtilcL  leave  the  eqaetims  Invariant.  (Note  the  requlrasent  p 

oy  the  presence  of  en  exterior  medl™.)  Ctie  easily  finds 


p  loosed 


Introducing 


e 

=  or  .  1 

6 

=  2(0-1)  . 

s  - 

?/t" 

— ♦ 

V  ^ 

P  = 

p/t^'-i) 
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the  equsfcicms  beccme 


'  V  +  (V  «  Qfs)  *  7p  =0 

pC(cif‘l  )V  *  (V-ofi)  .  #]  =  . 

2(«-l)  C;Pp"®  +  (V-as)  .  7(Pp^)3  ^  0 

All  the  quantities  are  now  functi<Mis  of  s  alme,  and  7  now  applies  to  s. 

For  siBiplieity  (tractahilityl )  one  now  searches  for  a  solution  ^ 
the  form 


V  =  A  +  r  s 

where  A  is  a  ccmstant  vector  and  F  a  e^stant  dyadic  (ii»trix),  Then 


7  •  V  =  trace  F  =  a,  say 


Let 


and 


s  ss  r  -  ffi 

U  =  V  .  a  1 

=  A  +  S  s 

the  equations  then  have  the  for® 


(A  +  Ss)  •  7p  +  ap  =0 

(r+(o-i)i)S  +  F(r-i)s  +  i  7  p  =  0 

p 

(A  +  Hs)  •  7(pp'’^)  +  2(cr-l)p  s~'^‘  -  0 

The  first  and  last  eqiiaticms  may  he  solved  hy  the  method  of  characteristics. 
Placing 


f  -  log  p/p^ 


s) 


(me  has 
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^  ^  1 

"  a 


U 


with 


s  =  s 


P  -  p„  ,  V  =  v„  ,  A  .  V.  -  rs.  . 


for  the  initial  conditiais.  Hence 


.  t  = 


U=V^-rs  +ss  =  e 
o  o 


a 


(Vg  -  • 


Similarly,  putting  x  =  p  "the  third  equation  gives 

U  =  e  (V  -  as  ) 

'  o  o 


whence 

X  -  (v  +  ♦ 

p 

I  =  (i)  “ 

p  'p  ' 

O 

=  (v  +  §i®lil)p  Yi.  . 

a  p 

Inserting  these  quantities  into  the  second  equation  OTie  finds  (eventually] ) 

e  )  1^ 


(1  .  p  “  )  +  oe-^  (1  .  e  ®  ^ 


u  . 


(r  +  (c^-i)i)  (v^  _  r  S^)  +  u  •  (r  -  1)  r  S  -  (ay  >  2(ff-i))  i  ==  o 

(A  +  £  s)  •  (r  +  (cr-l)l)  A  +  (A+Ss)  (F  -l)rs 

y  +  ,  1 


(ay  +  2(Qf-l))  (^) 


a 


o  "o 

-  A  •  (r  +  {ff-i)i)  A  +  A  [e  (r  +  (q!-i)i)  h  +  (r  -  1)  f]  s 

+  1  [s  (r  -  1)  r]  s 
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Since  H  is  a  3  dimensional  mtrix  (dyadic)  all  functicxis  of  S  n^y  be 

2 

written  as  the  sum  of  3  terms  ...  l+..,H+...Sj  the  coefficients 

—  f  %  /a  )c 

above  will  involve  the  proper  values  c.  of  and  the  functions  e  i^ 

^  A 

and  possibly  their  derivatives  up  to  seccsid  ordei'.  Since  p  =  %  the 

coefficients  of  terms  involving  derivatives  must  vanish.  Since  caily  one  term 
appears  on  the  right  of  the  equatioi,  the  bilinear  term  s  must  vanish  on 
the  left.  Hence  one  must  have 


H  (r  -  i)  r  =  0 


i.e.,  (r  -  oi)  (r  -  1)  r  =  0 

To  avoid  inessential  canplicaticais  for  the  present  let  7  =  F,  i.e.  F  is 
supposed  synanetric,  it^lying  V  X  D  =  0  a  reasonable  assumpticsi  for  many 
appli cat i  cns .  Thus 

(F  -  ol)  (F  ..  1)  F  =  0 


Hence  a  and/or  1  and/or  0  must  be  proper  values  of  F.  The  case  F 
corresponds  to  uniform  motion,  cuad  is  uninteresting.  The  case 


.  I 


F  = 


0 


corresponds  to  a  plane  situaticn,  the  case 

ill 


F  = 


to  a  cylindrical  (line  source)  configuration,  and  the  case  F  =  1  to  a 
spherically  syasnetric  a  point  cylindrically  symnetric  sitmtioi.  The 
presence  of  the  vector  A  (which  is  generally  independent  of  F)  allows  a 
superposition  of  uniform  and  expanding  (contracting)  motions.  The 
paiameter  or  maxst  be  determined  by  the  exact  shock  ccxiditions,  by  the 
properties  of  the  perturbation  soluticxi,  or  hy  the  integred  relatiais. 

The  specific  solutions  for  these  cases  are  given  hy  Raizer, 
and  their  ccmblnaticsi  with  conservation  laws  discussed  there. 


10 


III.  ASYMPTOTIC  BEdAVIOR;  PERTURBATION  S^UATIONS 


If  q)  is  a  vector  representing  the  dependent  varJables,  the  hydro- 
dynamic  equations  have  the  form 


at 

where  X  is  a  differential  operator  linear  in  the  derivatives.  A  similarity 
transformation  (of  the  type  considered  above)  implies  the  existence  of  a 
transformation  operator  S  such  that 

i  =  S"^  cp 


is  independent  of  t.  If 


then 


If  =  <p  -  S  § 


H  =  X  (cp)  -  X  (S§) 

=  X  (SI  -  Y)  -  X  (SI) 

^  X'  (Si)  ¥ 

plus  higher  order  terms.  These  first  order  perturbation  equations  have 
the  formal  solution 


t 

f  X'  (Si)  dt 

Y(t)=e‘°  f(o)  . 

Thus,  if  the  set  of  equations  admits  such  a  similarity  transformation  S, 
it  is  conceivable  that  a  similarity  solution  i  can  be  found  such  that 
¥(t)  0  as  t  -•  even  though  ¥(0)  ^  0  may  be  large. 

Even  if  the  equations  do  not  admit  an  exact  similarity  solution 
(e.g.,  due  to  equation  of  state  difficulties),  it  is  possible  that 
approximations  to  the  equations  do  so,  and  that  the  remaining  terms  may 
be  treated  as  an  inhcsnogeneous  perturbation,  which  under  certain  conditions 
-*  0  as  t  -  «>  in  analogy  to  the  right  hand  side  of  ¥(t)  above. 


To  conclude,  it  is  of  interest  to  write  down  the  perturbation 
equp lions  for  the  type  of  similarity  transfonnation  given  above.  Using 
the  notation  above,  one  writes 

p(t»r)  -  p(t,s)  =  R(s)  +  G(t,s) 

v(t,r)  -  v(t,s)  =  V(s)  +  w(t,s) 

P(“>r)  -  p(t,s)  -  p(3)  +  q(t,s) 

where  R(s)  has  now  been  written  in  place  of  p(e)-  The  perturbation 
equations  have  the  form 


where  H  =  {a,  w,  q) . 
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